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We study the motion of free magnetic polarons in a paramagnetic background of fluctuating local 
moments. The polaron can tunnel only to nearby regions of local moments when these fluctuate into 
alignment. We propose this fluctuation induced hopping as a new transport mechanism for the spin 
polaron. We calculate the diffusion constant for fluctuation induced hopping from the rate at which 
local moments fluctuate into alignment. The electrical resistivity is then obtained via the Einstein 
relation. We suggest that the proposed transport mechanism is relevant in the high temperature 
phase of the Mn pyrochlore colossal magneto resistance compounds and EuB6. 
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Recently Free Magnetic Polarons (FMP) have received 
renewed attention. They were proposed to explain the 
Colossal Magnetoresistajiae JjCMR) in the Manganese 
Pyrochlore compoundstrHuQ'ElEl and they have been stud- 
ied in the context of the double exchange,Hiodel and the 
Manganese Perovskite CMR compou»d«ErErcl. Moreover, 
Raman scattering data has suggestedt3 : til that they exist 
in EuBg. Previous theoretical studies have focussed on 
the static properties of the FMP. Here we focus on the 
dynamic aspect, propose a new transport mechanism for 
a FMP and calculate the resulting resistivity. 

A magnetic polaron is a composite object consisting 
of a localised charge carrier and the alignment it induces 
in a background of local moments. Localisation can oc- 
cur for two different reasons: the carrier can be trapped 
by an impurity atom and then induce a magnetisation 
in the region where it is localised. The resulting par- 
ticle is called a "Bound Magnetic Polaron" (BMP). It 
is well documented experimentally, for example in di- 
lute magnetic semiconductors like Cdi-^Mn-rSeEj, and in 
rare earth clpicogenideaiix It has been studied in depth 
theoreticallyE-1 A BMP is not free to roam through the 
sample since it is bound to its impurity. Only activated 
transport is possible: when the BMP is "ionised" the 
carrier is free to move until it is trapped by the next 
impurity. 

However, for large enough coupling to the local mo- 
ments the cafirifrr, can self-trap without the need for 



an impurityOEjO, forming a Free Magnetic Polaron 
(FMP). Due to the coupling the carrier acts as a mag- 
netic field on the local moments. The strength of this 
field varies in space as the probability density of the car- 
rier: the more localised the carrier the stronger the field 
and the larger the energy gain resulting from aligning the 
local moments. The region of aligned moments there- 
fore acts as a potential well that localises the electron 
and a FMP is formed. The balance between the gain 
in magnetic energy from induced alignment and the loss 
in kinetic energy because of localisation determines the 
polaron size. The existence of an FMP has not been con- 
firmed experimentally, but it has been suggested in the 



E 



E v 



\ ^r s ^ 






\ Level crossing 
and hop 




Thop 




TA 


TX log(t') 



FIG. 1: Time evolution of electron state. Black curved 
lines and arrows in the upper part represent the electron den- 
sity and spin. A carrier is put on the "blue" level. A hop 
occurs when the "red" and "blue" levels cross. The density of 
states g(E) is shown on the LHS of the lower part. 

Mn-pyrochloresBB! and in EuBg as 

The mechanism of transport by FMP is in doubt. The 
conventional view is that transport is necessarily acti- 
vated, as for a BMP. Here we present an alternative view- 
point: we propose that, unlike the BMP, the FMP can 
move between nearby sites without thermal assistance. 
We consider a fluctuating, paramagnetic background of 
local moments. A neighbouring region of local moments 
can fluctuate into the same alignment as the polaronic 
moments. At that moment, the carrier can tunnel to the 
newly aligned region without needing to overcome an en- 
ergy barrier. The tunneling process is fast compared to 
the spin fluctuations. After the tunneling process the car- 
rier and the alignment have moved so that the complete 
FMP has hopped to the new location. The entire time 
evolution of the polaron formation and hopping process 
is illustrated in Figure [jj We call this transport mecha- 
nism "Fluctuation Induced Hopping" (FIH). It does not 
involve an activated process. We calculate the resistivity 
FIH gives rise to. 

In the next section we present a model Hamiltonian 
that provides the frame of reference for our work. Wc 
describe the static properties of the FMP and the band 
states in Sec. El postponing a justification until Ap- 
pendices |X] and In Sec. Ill we calculate the electrical 
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resistivity for polaron hopping. We determine the rate at 
which nearby regions of local moments accidentally align 
themselves due to fluctuations. Since the FMP tunnels to 
these regions, this rate determines the diffusion constant 
and hence the electrical resistivity. 



I. MODEL HAMILTONIAN 

We consider a low density electron gas that is coupled 
ferromagnetically to a background of local moments. The 
local moments are themselves coupled ferromagnetically. 
The following Hamiltonian describes this systemM: 

H = t^TciaCjo ~ J'^2&i- Si- J ^2 Si- Sj. (1) 

(i,j)° » (id) 

Here i denotes the lattice site, c\ a creates a conduction 
electron, S% is the local moment on site i and &i is the 
conduction electron's spin. denotes a summation 

over nearest neighbours. The first term of the Hamilto- 
nian in Eqn. [l] is the kinetic energy of the carriers, the 
second term couples the carriers to the local moments 
on which they reside. The third term couples the local 
moments ferromagnetically. This term can be due to, for 
example, superexchange. We have an s-d Hamiltonian 
with an additional Hcisenberg term. We consider the 
strong coupling regime in which 



J' >t ~ 10J~ O.leV. 



(2) 



In our calculations we use the values J' = 5i and J = 
O.OleV which are in agreement with experimental values 
in the relevant materialsErE3. It should be noted that the 
magnetic transition is not driven by the s-d part of the 
Hamiltonian, but only by the superexchange because of 
the low carrier density. 



II. POLARON AND BAND STATE 

Here we present the wave function we use to study the 
transport properties of the FMP. Our variational calcu- 
lation in Appendix [A| shows that the FMP is small in 
the strong coupling regime: the carrier occupies approx- 
imately two lattice sites (see Figure ||) . We therefore use 
the following wave function to describe the polaron: 
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where Si are the polaronic local moments and m = 
So + Si is the magnetisation of the FMP. It describes a 
carrier localised on two lattice sites, with its spin quan- 
tised and pointing "up" along the direction of rh to min- 
imise the s-d energy. The s-d term in the Hamiltonian 
therefore reduces to —J' am, which is lower for more 
aligned local moments. This means that the carrier in- 
troduces an additional coupling between the polaronic 



moments. The value of the effective coupling constant 
can be obtained by expanding the expression for m up 
to first order in So • Si — S 2 for nearly aligned polaronic 
moments. We obtain 



Jeff = 



J' 



v/l + 8S(2S+l) 



2 J. 



(4) 



Since T <C J' the moments are nearly aligned, so the 
polaron energy is 
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In addition to the polaron state with its induced mag- 
netisation there are many more possible states for the car- 
rier in which it does not align any local moments. These 
are states in the narrowed band described in Ref. [l^. 
Since the background fluctuates, these "band states" per- 
sist at a given location only for a time comparable to 
the timescale of these fluctuations, which we denote rg. 
Nevertheless, the FMP would be unstable if a significant 
number of lower energy band states existed, since the 
carrier could then tunnel to them and gain energy. 

To check whether the FMP is stable we need to esti- 
mate the position of the band edge. If the polaron level 
lies below the band edge, band states with a lower en- 
ergy are exceedingly rare, and can be neglected. If on 
the other hand the band edge lies below the energy of 
the FMP, the latter is unstable. We determine the po- 
sition of the band edge as the lowest energy of a typical 
band state. We use a variational approach to calculate 
this energy (details are given in Appendix |^). The size 
of the band state is determined by the balance of the ki- 
netic energy cost of localisation and the gain from the s-d 
term. The latter is very small: the carrier aligns its spin 
with the total magnetisation of the region it is localised 
in. This magnetisation is due to statistical fluctuations 
in the paramagnetic regime, and consequently very small. 
We therefore expect the kinetic term to dominate and the 
band state to be very large. This is confirmed by our cal- 
culation. We show that typical band states are extended 
over a region of roughly 10 6 lattice sites. Their total en- 
ergy is very close to — 6t (See Equation (|B^)). Therefore 
the polaron level lies below the band edge in the regime 
we consider (J' > t), and the FMP is stable. This con- 
cludes our discussion of the static properties of the FMP 
and the band state. In the next section we determine 
their dynamic properties. 



III. FLUCTUATION INDUCED HOPPING 

In this section we consider the FIH mechanism in de- 
tail and calculate the hopping rate of the polaron and the 
resulting electrical resistivity. Let us examine the time 
evolution of a single carrier that is injected onto a lat- 
tice site in an empty system. The local moments in its 
vicinity cannot respond immediately to the carrier's pres- 
ence, but react on a timescale ta- On timescales smaller 
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than ta, the background appears static. A completely 
static background would cause Anderson localisationEa 
and trap the carrier in a state in the tail of the band 
below the mobility edge. Although localised, the car- 
rier has not yet induced any alignment between local 
moments in its vicinity; it is in one of the band states 
described above. Only when the moments have aligned 
themselves is the energy of the state greatly reduced. 
This means that there is a large energy barrier that pre- 
vents the carrier from making thermal hops out of the 
region of alignment. However, the energy of band states 
fluctuates because the local moments fluctuate. It can 
therefore cross the polaron level. At such a crossing the 
electron can tunnel to this level, since there is no more en- 
ergy barrier to overcome. We call this tunnelling process 
Fluctuation Induced Hopping (FIH). After the electron 
tunnelled, the entire FMP has moved: carrier and align- 
ment have jumped to another site. We summarise the 
entire time evolution in Fig. [l]. 

The occurrences of hops are uncorrelated in time and 
space since the background of local moments is param- 
agnetic. FIH is therefore a Markoff process and the FMP 
executes a random walk. In an ensemble of realisations of 
the polaron and the background, polarons in different re- 
alisations follow different paths. The Probability Density 
Function (PDF) of the polaron, defined as the fraction 
of realisations in an ensemble that have the FMP at a 
specified time at a specified location, obeys the diffusion 
equation.Ej The diffusion constant in this equation char- 
acterises the polaron transport in the long time limit. 

For a random walk in three dimensions that consists 
of hops of I lattice Constantsa-occurring with a frequency 
lui, the diffusion constant isnd 
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where a is the lattice constant. The resistivity from po- 
laron transport is then obtained from the diffusion con- 
stant by means of the Einstein formula: 



P = («e/i) 1 = 
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where /i is the mobility, and n the number density of 
polarons. 



A. Rate of level crossings 

We calculate the rate at which band state levels cross 
the polaron level. Such a crossing occurs when the band 
state energy fluctuates so much that it lies at or below 
the polaron level. The crossing rate depends on the size 
of the band state: the energy gap between the FMP and 
band states depends on their size. In addition to this, we 
will see that the energy of large band states fluctuates 
less. First we calculate the crossing rate for small band 
states. 



The energy difference between the polaron and a small 
band state is only due to the difference in exchange en- 
ergy since their kinetic energies are identical. The time 
at which the RMS deviation of the band state energy be- 
comes as large as the energy gap between the two levels 
therefore determines the crossing rate. 



({So ■ Si)j eg — (So ■ Si)j) 
([So(tk) ■ Si(tx) - So(0) • Si(0))j] 2 ) 



(8) 



where angular brackets denote ensemble averaging. We 
neglect the variance of the polaronic exchange energy 
since it is very small due to the large effective coupling 
constant of the polaronic moments (See Appendix ^|). 
We also neglect the exchange energy of the band state 
compared to the polaron exchange energy: 



[So(tx) ■ Si(t x )][S ■ Si]) j 
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For simplicity we use an interpolation between t = and 
t = oo instead of the exact analytical form of the four- 
point correlator in Eqn. ^. At time t = the correlator 
reduces to ((S ■ Si) 2 ). A pair of spins at t = is com- 
pletely uncorrelated with itself at t = oo, so that in this 
limit the correlator reduces to (So- Si) 2 . We can therefore 
interpolate as follows: 



([S (t) ■ Si(t)][So(0) ■ Si(0)]) 
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TS 



((So ■ Si) 2 ) - (So ■ Si)' 
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where f(x) varies smoothly from at x = to 1 at 
x = oo. Moreover, we have written /(t/rg) since the 
four-point correlator varies on the same timescale as the 
fluctuations of the background. This allows us to rewrite 
Eqn. [)] as: 
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((So • Si) 2 ) j ~ (So ■ Si 



(11) 



Since f(x) is a smooth function varying between and 1 
that changes mostly near x = 1, f(l) » 1/2 and f'(l) ~ 
1. We expand /(t/rg) up to first order about t = rs, 
which yields 



(So -S 



((S -Si) 2 )j- (So -Si) 2 



(12) 



This means that crossings occur on the time scale of the 
fluctuations of the Heisenberg magnet weighted by the 
different alignments of the polaron and the background 
spins. Since fcgT <C J' , the numerator of Eqn. 12 reduces 
to (S(S + 1)) 2 . 

For k^T > O(J) the two terms in the denominator of 
Eqn. [l^ have very similar temperature dependences, so 
that 



T X = At S 



(SQg + l) 2 
(So -Si) 2 



(13) 
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where A ~ 0(1). We use the expression given in Ref. |2l] 
in the denominator and ts = h\J (3/ J.E3 to obtain 

h 



tx 



(14) 



Tx increases with temperature. With increasing temper- 
ature level crossings become more rare. The decrease in 
the timescale of the spin fluctuations is more than offset 
by the increase of the average misalignment of the local 
moments. 

For fceT ^ J the temperature dependences differ: 
((Sq ■ Si) 2 } tends to S 2 /3 since the spins are com- 
pletely uncorrelated, whereas (So • Si)j vanishes as 
Moreover the times cale of the fluctuations is different E3 
Tg = h/y / S(S + 1) J. Hence tx tends to a constant in 
this limit: 

r x = / 3h . (15) 

Jy/S{S+l) 

The reason is that the local moments are completely dis- 
ordered in this regime; an increase in temperature does 
not cause an increase in disorder. The crossover be- 
tween the two regimes occurs at a temperature of about 
JS(S+l)/k B . 

Equations (jl4|) and ([H]) constitute the principal re- 
sult of this section. However, before proceeding to the 
estimate of the diffusion constant we should check that 
fluctuations of other band states - in particular those in- 
volving rearrangements of many spins - do not change 
our conditions. We calculate the crossing rate of large 
band states with the polaron level in the same way as 
before (see Appendix ^| for details). For a level cross- 
ing with a large band state, many local moments need to 
fluctuate into alignment simultaneously. This is a very 
unlikely event, and one expects the crossing rate to be 
accordingly small. This expectation is borne out by our 
calculation. In Appendix |C| we show that level crossings 
with large band states can be neglected safely. 



B. Diffusion Constant 

When a crossing occurs it is possible, but not neces- 
sary, for the FMP to hop. The probability, Vi, of a hop 
of length I at a level crossing depends on the overlap be- 
tween polaron and band state wave functions and on the 
rate at which the levels cross. The frequency at which 
hops of length I occur, uji, is then 



(16) 

since the number of small band states a distance I away 
from the FMP is roughly proportional to I 2 . The tun- 
nelling probability from the polaron state, |P), to the 
band state, \B) with energies Ep and Ep, is given by the 
Landau-Zener formulaE-JCJ 



1 — exp 



2tt \(P\H\B) 



^(Ey 



i at 



E B ) 



(17) 



where H is the Hamiltonian for the particle. Here we 
have \-^(Ep — E B )\ ~ J'S/ts, since the energy difference 
is due to the initially unaligned local moments in the 
band state. The time in which this difference between 
the levels disappears is rg . The spatial extent of the wave 
functions of the polaron and the small band state limits 
the hopping range to one lattice constant. The overlap 
between neighbouring small polaron states given is by: 



(A\H\B) = ±(0\(co + Cl )H(cl + ci)\0) 



J' 



(18) 



where the sites and 1 are nearest neighbours and 1 
and 2 are nearest neighbours. Therefore the hopping 
probability at a level crossing between two neighbouring 
energy levels is given by: 
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(t + J'^f2S(2S 
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since 1 <C J'/ J. This means also that the probability 
is largely independent of temperature. We will therefore 
take the hopping probability at a level crossing between 
two neighbouring levels to be 1. The diffusion constant 
and resistivity for FIH are therefore 
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(20) 



This is our main result. We plot the resistivity versus 
temperature in Figure |^ interpolating between the high 
and the extremely high temperature regime. Firstly, 
the FMP can only hop to neighbouring sites when a 
favourable statistical fluctuation aligns the local mo- 
ments. These fluctuations are statistically likely in the 
sense that they can be estimated from the RMS devia- 
tion of the spin fluctuations. Occurrences of alignment 
far from the polaron do not lead to hopping. Secondly, 
while small polarons typically hop by this process, large 
ones cannot. The required spin fluctuation into the cor- 
rect configuration is statistically very rare. Thirdly, the 
timescale of the spin fluctuation is slow enough that the 
FMP hops with probability 1 once the requisite con- 
figuration is obtained. Fourthly, the diffusion constant 
decreases and the resistivity increases as a function of 
temperature. This reflects the increasing time intervals 
between the level crossings for higher temperatures and 
the relation between the resistivity and the diffusion con- 
stant. 

The resistivity we obtained is "metallic" : it increases 
with temperature, even though we are not considering 
a metallic system at all. It is interesting to compare 
our result to the resistivity of a very dirty metal, where 
fcpAf w 1, Af being the mean free path for the carriers. 
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FIG. 2: Temperature dependence of the resistivity. 

In such a material the Drude formula for the resistivity 
yields: 



It is clear that despite its temperature dependence the 
polaron hopping resistivity is far too large to be confused 
with scattering of metallic carriers. The mean free path 
in the dirty metal would need to approach one lattice 
constant for the resistivities to be comparable. At such 
short mean free paths the metallic picture of delocalised 
carriers breaks down. 



IV. CONCLUSIONS 

We have proposed a new transport mechanism for 
FMP in a fluctuating disordered background of local mo- 
ments. In our theory the FMP hops at the occurrence 
of a favourable fluctuation. The transport mechanism is 
therefore not activated, but gives rise to a "metallic" re- 
sistivity. Experimentally our theory can be checked by 
a measurement of the temperature dependent resistivity. 
Such measurements have been performed in two systems 
in which the presence of a FMP has been suggested: the 
Mn-pyrochlores and EuB*l 

In the Mn-pyrochloresoc3E3 the resistivity decreases 
with increasing temperature above T c . This is not in ac- 
cord with our predictions. There are several reasons for 
this. There arc different types of disorder that affect the 
resistivity, but were not taken into account in our work. 
In the In-doped compounds^!, there is a miscibility gap 
for dopings of 0.5 < x < 1.5. In this regime the bulk 
consists of two types of grains, each of a distinct phase 
with a different lattice parameter. Transport is domi- 
nated by processes associated with the grain boundaries. 



Phase separation is suspected to occur in the Sc-doped 
materials as well. A different type of disorder occurs 
in Bi-doped compounds: a Bi ion introduces a strong- 
scattering 6s vacancy on the Tl-sublattice. This could 
bind the polaron, making its hopping activated. Both 
the scattering centres and the grain boundaries make the 
predicted resistivity difficult to observe. 

EuBg is a much cleaner system, in which magnetic po- 
larons have been .observed by spin flip Raman scatter- 
ing experiments I 1 ™ 11 ! There is good qualitative agreement 
with experiment in the high temperature paramagnetic 
regime. The resistivity increases rapidly with tempera- 
ture up to about 150K and then more slowly; at tem- 
peratures above 200K the resistivity increases even more 
slowly.EScJ This agrees qualitatively with the crossover 
we predict. The crossover temperature we predict, 225K 
is only a little bit too large. The slowdown above 200K 
is probably due to other scattering mechanisms for the 
local moments, such as spin-orbit coupling and scatter- 
ing by carriers. There is no quantitative agreement as no 
T 5 / 2 law is observed. This is probably due to material 
specific complications that our theory does not take into 
account: EuBg exhibits two distinct magnetic transitions 
between phases with different magnetic anisotropics Ej 

There is good agreement as well with the 2-dimensional 
Monte Carlo simulation^. The simulation and our work 
agree qualitatively on the static characteristics of the po- 
laron, such as the temperature and J' dependence of its 
size and the core magnetisation. We also agree on the 
temperature window in which the FMP exists. We do 
not expect more than qualitative agreement given the 
different parameter regimes that were explored: the sim- 
ulation uses a much weaker superexchange coupling. The 
discrepancy between the results for the binding energy in 
the critical regime can be understood. The simulations 
show a decrease in magnitude with temperature and we 
predict an increase. This is due to the breakdown of our 
high temperature approximation. This also explains why 
the simulation observes a larger polaron: we neglect the 
correlations between local moments except those induced 
by the presence of the carrier, whereas the simulation 
takes all correlations into account. 

The results for the dynamic characteristics of the FMP 
are different. The simulation shows a diffusion constant 
that is nearly independent of temperature. This result 
does not take into account the decrease of the timescale 
of the fluctuations: the diffusion constant is calculated 
as the average hopping per unit Monte Carlo time, de- 
fined as the number of rediagonalisations. To correct 
for this we divide the numerical diffusion constant by 
T§. This correction increases the discrepancy since the 
diffusion constant now decreases as a function of tem- 
perature. The discrepancy between the corrected Monte 
Carlo diffusion constant and our result indicates that the 
simulation looks at a slightly different form of polaron 
transport. In the simulation the polarons are larger than 
in our work. Larger polarons can move if a shell of neigh- 
bouring local moments aligns themselves. This motion 
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"through accretion" was not taken into account in our 
result, which was derived for a small FMP. Our theory 
and the Monte Carlo simulation do therefore agree qual- 
itatively where agreement is expected. 
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APPENDIX A: SMALL SPIN POLARONS 
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We determine the size and energy of the FMP in strong 
coupling regime by a variational calculation. We choose a 
trial wave function for the carrier and obtain the electron 
density at the local moments, which acts as an external 
field on the moments. The resulting alignment and de- 
crease in energy is calculated using Curie- Weiss theory. 
The expectation value of the energy of the trial wave 
function is then minimised with respect to the size of the 
wave function. 

A Gaussian-like function is used for the electronic part 
of the trial wave function for the FMP; in the notation 
of Eqn. |: 

l P > = E 6-^4, |0) ® |T) M <8> |m) , (Al) 

fi 

where the vectors r, are the positions of the lattice sites, 
measured from the centre of the polaron and Af ensures 
proper normalisation. |T)^ denotes the electrons spin, 
which is quantised and pointing "up" along the direc- 
tion of the average magnetisation induced by the carriers 
presence. |m) is the state vector of the polaronic local 
moments. The wings of this wave function take into ac- 
count that a trapped electron can nearly always make 
short excursions to a neighbouring non-polaronic local 
moment. This is possible since it has nearly always a 
spin component parallel to this moment. These excur- 
sions diminish the polaron energy insofar as they reduce 
the magnetisation of the core through a reduced electron 
density. 

The magnetisation of the background resulting from 
the presence of the carrier is obtained from Curie- Weiss 
theory: 

m(r) = ^B | [/3(Jm+^p ft )] ! (A2) 

where B is the Brillouin function and is the electron 
density at site n. Curie- Weiss theory neglects spatial 
correlations between the local moments that are not due 
to the effective field of the carrier and the induced mag- 
netisation itself. These spatial correlations are small in 



FIG. 3: Polaron Size vs. Temperature. The curves stop 
when the polaron energy reaches the edge of the band (esti- 
mated in Sec. ^). 



a paramagnet, since the coupling constant of the non- 
polaronic moments is much smaller than the effective 
field in the core of the FMP (J < J'). Curie-Weiss 
theory is therefore accurate in the core of the polaron. 
In the critical regime non-polaronic correlations become 
important. Our theory is not valid there. 

The expectation value of the energy of the trial wave 
function is minimised numerically with respect to the 
size, A, of the polaron. The results are shown below in 
Fig. |. We see that the FMP is small: its size is of the 
order of a lattice constant over a wide range of temper- 
atures for realistic values of the parameters. The calcu- 
lated magnetisation is always close to saturation in the 
core of the polaron, but decreases slightly as a function of 
temperature. For higher temperatures the FMP shrinks 
and its energy increases. This is because at constant size 
the energy gain from the s-d term decreases as the tem- 
perature is raised. To compensate the polaron shrinks, 
so that the effective field due to the carrier increases, 
which results in a magnetic energy gain. Hence, at higher 
temperatures, the point of minimum energy is shifted to- 
wards smaller sizes. This also means that a smaller s-d 
coupling increases the size. The increase of the polaron 
energy with temperature can be understood as follows: 
the kinetic energy is independent of temperature and a 
decreasing function of size. The magnetic energy curve 
shifts up as temperature is raised and is an increasing 
function of size. Therefore, the minimum value of the to- 
tal energy increases with temperature. These results are 
in close agreement with Ref. [l] which refers specifically 
to the pyrochlores. 

Our calculation is well-behaved at the ferromagnetic 
transition, even though this temperature lies outside its 
range of validity. As T c is approached from above, the 
polaron size tends to a value of approximately — 9At at 
T c . Neither the present calculation, nor Ref. [I] takes the 
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correlations of local moments outside the FMP into ac- 
count. The predicted size is therefore too small in the 
critical regime and the energy is overestimated. 



APPENDIX B: BAND STATES 

We check whether the FMP we discovered in the pre- 
vious section is bound by comparing its energy with the 
position of the band edge. We determine the position of 
the band edge as the lowest energy of a typical "band 
states" . The energy of a band state is the difference be- 
tween the energy of the system with an electron present in 
that state and the energy of the system in the absence of 
that electron. Again we use a variational approach with 
a trial wave function similar to Eqn. [Al[ with the caveat 
that local moments are not aligned. Only the kinetic and 
the s-d term of the Hamiltonian in Eqn. ^ contribute to 
the energy, the Heisenberg term is the same regardless of 
the presence of the carrier in a band state. The kinetic 
energy of a localised state is determined as before. The 
magnetic energy is estimated as follows. There is no net 
magnetisation in the paramagnet, so the s-d term in the 
Hamiltonian vanishes on average. In a finite region how- 
ever, there are statistical fluctuations that make it devi- 
ate from this average value, so that there is a non-zero 
magnetisation in this region, rha uct = J2 i Pn^fi- The 
carrier's spin is quantised along the direction of mfl uct , 
pointing "up" . The average value of the s-d term in the 
Hamiltonian is then 



J' 



(Bl) 



The above sum is split up in one where i — j and one 
where i ^ j. We then take the thermal average of the 
Taylor expansion of the square root about the i = j term. 
The sum that contains term with i ^ j is at least of order 
f3J since it contains correlations between different local 
moments. The s-d energy of the band state is then given 
by 



■0(fiJ)]. 



(B2) 



The energy of the band state is minimised numerically 
with respect to the size. We obtain a size of the order of 
10 6 local moments that decreases weakly with increasing 
J'. We find that the s-d energy is completely negligible 
and that the band state is so large that its energy is 
very close to — 6|i|. We expand the band state energy to 
second order in 1/2 for large / and obtain: 



£Band = -6|i|(l-i)-^ 



1 



(V21) 3 



(B3) 



where I is the spatial extent of the band state. The po- 
laron level lies therefore below the band edge and the 
FMP is well bound. 



APPENDIX C: CROSSING RATE WITH LARGE 
BAND STATES 

We use the same method to calculate the crossing rate 
for large band states. However, instead of two spins fluc- 
tuating into nearly exact alignment many spins need to 
collectively fluctuate into a more aligned configuration. 
The time between two crossings is again estimated from 
the time it takes for the RMS deviation of the band state 
energy to become as large as the energy difference be- 
tween the band level and the polaron level: 

(CI) 



{E n [tx(N)\ - E N (0)Y ) = (E P - E N ) 



Here tx{N) is the average time between crossings of the 
polaron level and a particular band state of N local mo- 
ments with energy En. In Appendix [B| we calculate en- 
ergy of the band state using a variational Ansatz with an 
accurate Gaussian trial wave function to show that the 
FMP is bound. Here, however, we determine the cross- 
ing rate and we will see that a simple model is sufficient. 
We assume that the carrier is localised uniformly with- 
out inducing extra alignment. Its spin quantised along 
the direction of the sum of all the local moments in the 
band state, m, pointing "up" . The kinetic energy of the 
band state is approximately — 6|t|(l — tt 2 /N 2 ^) and its 
exchange energy vanishes on average. The kinetic en- 
ergy is c onstant in time, so it cancels out in the LHS of 
Eqn. CI, leaving only the RMS deviation of the exchange 
energy. 

Now we use the expression for the 2-site polaron of 
Eqn. | on the RHS of Eqn. |c|. The band state s-d 
energy can be expressed as am(t) /N 2 and hence Eqn. CI 
becomes: 



N 



{ [SfM ■ s rl (t)][s fk (o) • £ n (o)]}* 

\ i,j,k,l 



t=rx (jV) 



t=0 



2N l 



(5 



t 



N 2/3 I J, 



S 



-i 2 



) 11 9 



(C2) 



The sum on the LHS of this equation contains N 2 terms 
of the form S 4 and 2N terms of the form S 2 Y^Lj ' ^ 
where the spin operators are evaluated at equal times. 
There are also terms where i ^ j and k ^ I. The square 
root is expanded about the term of order N: 

N 

{ [SrM ■ sv- (*)][&» (0) • &,(0)]}* = 



N 



N 



NS(S+l) + J2Sr k -S Pl 

i¥=3 



Thermal averaging both sides of this equation results in 
many two- and four-point spin correlators. We only re- 
tain the correlators of lowest (quadratic) order in the 
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small parameter /3J, thereby considering only nearest 
neighbour interactions. Moreover, the first and second 
terms on the RHS of Eqn. cancel in Eqn. C2 since 



they do not depend on time. Hence the condition for a 
level crossing reduces to 



(S ■ Si) 2 - {5 [r x (iV)] • Si[tx(A0]}[So ■ Si 



2S(S+1)N 



2 r 



(5- 



Jt[ 

^2/3 > J, 



(C4) 



The four-point correlator is treated as in Eqn. |l0| a nd we 
introduce the numerical constant A as in Eqn. [13c 



/(rxW) = AN- 



(5 



25(5+1) 
3(5 • 5i) 2 

N 2/3 > J/ 



(C5) 



f(t) could be expanded about t = rg to obtain an 
explicit expression for rx(N). However, it is clear that 
large band states rarely cross the polaron level. f(t) on 



the LHS is bounded from above by 1 and the RHS is 
proportional to iV 2 , and so only small band states sat- 
isfy Eqn. |C5|. Since rjv grows with N, the size of the 
largest crossing band state follows from Eqn. C4 in the 
limit tx(N) — > oo. In this limit f(t) — 1 so that the 
largest crossing band state needs N ~ (So ■ Si) 2 <C 1. 
This result shows that large band states do not cross the 
polaron level according to Gaussian statistics. The phys- 
ical reason is that the contribution to the s-d energy of a 
single pair of local moments is weighted by the electron 
density. For large states this density is low, so that align- 
ment of a single pair of local moments does not lower the 
energy significantly. Thus, a crossing requires all the lo- 
cal moments to fluctuate into nearly perfect alignment. 
Of course, large clusters of ferromagnetically aligned lo- 
cal moments do exist, and do cross the polaron level, but 
they lie in the far tail of the band and are far more rare 
than a Gaussian approximation to the density of states 
would predict. These fluctuations are thus negligible. We 
have therefore shown that our expressions for rx in Eqn 
( |l4| ) and ( |l5| ) give the correct hopping rate for the small 
FMP. 
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